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Abstract 

We investigate the evolution of the global string network in the radiation 
dominated universe by use of numerical simulations in 3+1 dimensions. We 
find that the global string network settles down to the scaling regime where 
the energy density of global strings, p s , is given by p s = £/u/i 2 with p the 
string tension per unit length and the scaling parameter, £ ~ (0.9 — 1.3), irre- 
spective of the cosmic time. We also find that the loop distribution function 
can be fitted with that predicted by the so-called one scale model. Con- 
cretely, the number density, ni(t), of the loop with the length, is given by 
ni(t) = v/[t 3 / 2 (l + Kt) 5 / 2 ] where v ~ 0.0865 and k is related with the Nambu- 
Goldstone(NG) boson radiation power from global strings, P, as P = up with 
k ~ 0.535. Therefore, the loop production function also scales and the typical 
scale of produced loops is nearly the horizon distance. Thus, the evolution 
of the global string network in the radiation dominated universe can be well 
described by the one scale model in contrast with that of the local string 
network. 



I. INTRODUCTION 



Cosmic strings could be formed in consequence of the cosmological phase transition at 
the very early universe [0. They are divided into global and local(gauged) strings according 
to the property of the broken symmetry. Though they have similar properties, such as they 
are line objects with the false vacuum energy, intercommute at crossing [fj,[3| and so on, 
there are also many differences between them. While a local string has two cores comprised 
of a magnetic flux core and a scalar field core, a global string has only a scalar field core. 
For local strings, the gradient energy of the scalar field can be canceled out by the gauge 
field far from the core so that the core is well localized and the vacuum energy of the core 
is dominant. Therefore, the Nambu-Goto action is adequate to follow the evolution of the 
local string network except at crossing. On the other hand, there is no gauge field in the 
global string model so that the total energy of global strings is dominated by not the vacuum 
energy of the core but the gradient energy of the scalar field, namely, the NG boson field. So, 
we need consider not only the core but also the associated NG boson field and the coupling 
between them in contrast with the case of local strings. Thus we must use the Kalb-Ramond 
action instead of the Nambu-Goto action as an effective action M. 

The dynamics of the local string network has been examined by use of the Nambu- 
Goto action. Kibble || first proposed the so-called one scale model, where the behavior of 
the system can be characterized by only a parameter, namely, the scale length, L, and an 
unknown loop-production function. He showed that either the local string network goes into 
the scaling regime where the scale length, L, grows with the horizon distance or, the 
universe becomes string dominated. In the scaling regime, infinite strings intercommute to 
produce closed loops so that at any time the number of strings stretching across the horizon 
distance within each horizon volume is almost constant and produced loops decay through 
radiating gravitational waves |7J]. Bennett developed the Kibble's one scale model and 
showed that unless most of produced loops self-intersect and fragment into smaller loops 
with the typical length smaller than L, the reconnection rate is large enough to prevent 
scaling. Mitchell and Turok || studied the statistical mechanics of the string network in the 
flat spacetime and found that the equilibrium distribution of the string network is dominated 
by the smallest loops allowed, which suggested that strings tend to break into very small 
pieces in the expanding universe. Albrecht and Turok []TU| modeled the network as the hot 
body radiation where loops are radiated from the long strings and showed that the scaling 
solution is inevitable. However, the application of the flat spacetime statistical mechanics to 
the string dynamics in the expanding universe is not necessarily justified. Thus, numerical 
simulations are unavoidable to decide whether the local string network settles down to the 
scaling regime or not. 

Albrecht and Turok JTTJ] gave the first numerical investigation for the scaling property. 
Later three groups flIT?|P^| , |T3[ improved numerical codes and all groups concluded that the 
large scale behavior of the local string network goes into the scaling regime with the scaling 
parameter f equal to (50 ±25) |TD[, (13 ± 2.5) 0, and (16 ±4) fT3[ in the radiation- 
dominated era. Though Albrecht and Turok plj found that the loop distribution function 
also scales, the higher resolution simulations performed by the other two groups |l2|Jl3|] 
showed that long strings have significant small structures and loops are typically produced 
at scales much smaller than the horizon distance, which is close to the cut-off scale. In 
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response, Austin et al. |14| proposed the three scale model with £, the step length £, and 
the scale ( describing the small scale structure. They found that £ and £ grow with the 
horizon distance but ( begin growing only when the gravitational back reaction becomes in 
effect with C/£ ~ 10~ 4 . However, Vincent et al [HJ recently claimed from the flat spacetime 
simulations that instead of loop production due to intercommutation, long strings directly 
emit massive particles so that the dominant energy loss mechanism of local strings is not 
gravitational radiation but particle production. 

Thus, though there is a consensus that the large scale structure of the local string network 
obeys the scaling solution, the loop production and the dominant energy loss mechanism are 
now in dispute. This is mainly because inclusion of gravitational radiation to the numerical 
simulations is impossible and gravitational radiation is so weak that kinks live for a long 
time. 

On the other hand, the evolution of the global string network has been less studied. 
Its manifestation, however, is essential to estimation of abundances of relic axions radiated 



from axionic strings [pLq— |1S|| , which may be the cold dark matter. Also, scaling property is 
indispensable for the scenario where global strings become the seed of the structure formation 
of the universe and produce the cosmic microwave background anisotropy |19|]. Thus, it is 
important to clarify the dynamics of the global string network, in particular, whether it 
enters the scaling regime like the local string network. 

So far, the result for the local string network has been applied to the global string 
network because global strings also intercommute with the probability of the order unity 
@. But global strings have the associated Nambu-Goldstone bosons, which lead to long- 
range forces between strings and become the dominant energy loss mechanism of global 
strings. Therefore, instead of the Nambu-Goto action we have to use the Kalb-Ramond 
action [Q, which is comprised of three components, the Nambu-Goto action, the kinetic 
term of the associated Nambu-Goldstone fields, and the coupling term between them. But 
the Kalb-Ramond action has difficulty of logarithmic divergence due to self-energy of the 
string. Dabholkar and Quashnock |20| solved this difficulty by the similar prescription as 
given by Dirac in the electromagnetic system ||21|| , where divergence is renormalized into the 
electron mass. They gave the renormalized equation of motion for a global string comprised 
of the free part derived from the Nambu-Goto action with the damping term, which becomes 
negligible for a circular loop in the cosmological scale where \n(t/5) ~ (9(100). Then they 
concluded that the global string network can be well approximated by the motion of the 
Nambu-Goto action. However, as shown by Battye and Shellard |2^] though the calculation 
is done in the flat spacetime, the kinks on long string are substantially rounded due to the 
backreaction of NG boson radiation, which may significantly affect the small scale structure 
(if at all) of the global network system. Also, the above approach cannot include the long- 
range force between long strings, which may decrease the energy density of long strings. 
Thus, the examination of the dynamics of global strings by use of the Kalb-Ramond action 
is not yet complete. 

In the previous paper ||23|| , we manipulated the equation of motion of the complex scalar 
field instead of the Kalb-Ramond action. Then, we showed that the large scale behavior of 
the global string network goes into the scaling regime and £ for the global string network 
becomes 0(1), which is significantly smaller than that for the local string network. 

In this paper, we give a comprehensive analysis of the evolution of the global string 
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network based on the model adopted in [23|. We show the scaling of the global string 



network under more general situations and investigate the dependence of £ on boundary 
conditions and some parameters. Then, the loop distribution function is given in order to 
decide whether the small scale structure exists like the local string network. 

The paper is organized as follows: In the next section, we give the formulation of numer- 
ical simulations. In the section 3, we give the method of the identification of string segments 
and closed loops. Then, the scaling parameter £ and the loop distribution function are given. 
Finally, we discuss our results and give conclusions. 

II. FORMULATION OF NUMERICAL SIMULATIONS 

We consider the following Lagrangian density for a complex scalar field 

£[$] = g^ v {d^){d v ^ - V eS [$,T] , (1) 

where g^ v is identified with the Robertson- Walker metric and the effective potential V e s [$, T] 
is given by 

K ff [$, T] = ^($$t - v 2 ) 2 + ^T 2 $$t . ( 2 ) 




FIG. 1. One-loop finite temperature effective potential 14ff[<I>,T] of the complex scalar field. 

For T > T c = v^3?7, the potential \4g[$,T] has a minimum at $ = 0, and the U(l) 
symmetry is restored. On the other hand, new minima |$| m m = V\f^ ~ (T/T c ) 2 appear and 
the symmetry is broken for T < T c (Fig. [I]). In this case the phase transition is of second 
order. 

The equation of motion is given by 

+3H<l>(x) - t^V 2 $(x) = -KffM , (3) 
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where the prime represents the derivative d/d& and R(t) is the scale factor. The Hubble 
parameter H = R(t)/R(t) and the cosmic time t are given by 

2 _ _8v^^ _ J_ _ 

~ 3*^30^ ' 2tf ~ T 2 ' [ } 

where m v \ = 1.2 x 10 19 GeV is the Plank mass, is the total number of degrees of freedom for 
the relativistic particles, and radiation domination is assumed. We define the dimensionless 
parameter ( as 

^ rj (l67r 3 g*?7 2 ) ^ 

In our simulation, we take ( = 10, which corresponds to rj ~ (10 15 — 10 16 ) GeV with 
g* = 1000 and later investigate ( dependence on the result. The energy density at each 
lattice point is written as 

p(x) = &(x)&(x) + ^V$(x) • V&(x) + V eS [$,T] . (6) 

We take the initial time ti = t c /A and the final time tf = 75 1; = 18.75 t c , where t c is the 
epoch T = T c . Since the U(l) symmetry is restored at the initial time t = t i: we adopt as 
the initial condition the thermal equilibrium state with the mass squared, 



m 2 = 



d 2 V eS [\%T] 



dm 



(7) 

|$|=0 



which is the inverse curvature of the potential at the origin at t — ti. In the thermal 
equilibrium state, $ and <£> are Gaussian distributed with the correlation functions, 



P\*(x)&(y)\p) = I *\ , 1 coth ^ 2 o +m2 e* fc ^) 

V 1 ^ >U /equal-time J 2?T 3 o Ju 2 , ™2 2 



( 2 -) 3 2v ^ 



+ m? 



11 V ; Vy;IM /equal-time J (2tt) 3 2 2 



(8) 



(9) 



The functions $(£c) and $(y) are uncorrelated for x ^ y. We generate these fields for the 

initial condition in the momentum space, because the corresponding fields $(fe) and $(fc) are 
uncorrelated there. Then these fields are inverse Fourier transformed into the configuration 
spaces by the FFT formalism. 

Hereafter we measure the scalar field in units of t~ l , t and x in units of ti, and the energy 
density in units of t^ A . The equation of motion and the total energy density are given by 
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TABLE I. Five different sets of the simulations under the periodic boundary condition 





lattice number 


lattice spacing 
(unit = tiR(t)) 


C 


realization 


e 


(1) 

(2) 
(3) 
(4) 


128 3 
64 3 
256 3 
256 3 


2^/3/25 
4^/25 
V^/25 
2^/25 


10 
10 
10 
10 


300 
10 
10 
10 


0.90 ±0.06 
0.90 ±0.05 
0.99 ±0.09 
0.97 ± 0.04 




(7) 


128 3 


2^6/25 


5 


10 


0.88 ±0.07 


<E>(a;) 


+ j t ^) - ^v 2 $(x; 


) = f |$| 2 + ^ 


--£>, 

144 y 




(10) 


p(x) 


= $(x)&(x) + -V$( 


» ■ V&(x) + i f 


|$|2_H^ _ 

11 144/ 


c 2 

l s |<J>| 2 
^ 36t 1 1 ' 


(11) 



where A is set to unity for brevity. The scale factor R(t) is normalized as R(l) = 1. 

We perform the simulations in four different sets of lattice sizes and spacings (See TABLE 
|.): (1) 128 3 lattices with the physical lattice spacing Sx phys = 2V3tiR(t) /25. (2) 64 3 latt ices 
with 5x phys = 4:y/3tiR(t)/25. (3) 256 3 lattices with Sx phys = v / 3* i -R(t)/25. (4) 256 3 lattices 
with 5x phys = 2y/3tiR(t)/25. In all cases, the time step is taken as 5t = O.Oltj. In the case 
(1), the box size is nearly equal to the horizon volume (if -1 ) 3 and the lattice spacing to 
a typical width § ~ 1.0/(v / 2^) of a string at the final time tf. Furthermore, in order to 
investigate the dependence of £, we arrange the case with ( = 5, (7) 128 3 lattices with the 
physical lattice spacing Sx phys = 2y/6UR(t)/25. In this case we follow the time development 
of the system until the final time, tf = 150 tj = 37.5 t c , when the box size is nearly equal to 
the horizon volume (H -1 ) 3 and the lattice spacing to a typical width of a string. For each 
case, we simulate the system from 10 ((2), (3), (4), and (7)) or 300 ((1)) different thermal 
initial conditions. Since the simulation box is larger than the horizon volume even at the 
final time of the simulation, we adopt the periodic boundary condition. But, under the 
periodic boundary conditions, there exists no infinite string so that it is possible that string 
completely disappears in the simulation box. Therefore, in order to verify the dependence of 
boundary conditions, we also simulate the cases in TABLE |I| under the reflective boundary 
condition where V 2 $(x) on the boundary points disappears f\ 

Using the second order leap-frog method and the Crank-Nicholson scheme, the discretized 
equation of motion reads 



*Note that this condition is different from the open boundary condition where V<5(a;) on the 
boundary points disappears. Under the open boundary condition, the string feels attractive forces 
from the boundary so that the number of the strings tends to decrease than that in the real universe 
as under the periodic boundary condition. 
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TABLE II. Six different sets of the simulations under the reflective boundary condition 





lattice number 


lattice spacing 
(unit = UR(t)) 


>- 

c 


realization 




(1) 


128 3 


2^/25 


10 


10 


1.77 ±0.03 


(2) 


64 3 


4^/3/25 


10 


10 


1.57 ±0.04 


(3) 


256 3 


V3/25 


10 


10 


2.00 ± 0.05 


(4) 


256 3 


2^/25 


10 


10 


1.30 ±0.02 


(5) 


128 3 


4^/3/25 


10 


10 


1.18 ±0.03 


(6) 


256 3 


4^/25 


10 


10 


1.03 ±0.03 



l,ra+l/2 



1 , 3(5t 

4t L 



3 5 1\ ■ St 9 f , 

1 <£>• 1/9 ±-V 2 $.- -<5M $ 



C 2 

*' nl 36 1 



ill $2. 
144 | *.« 



V 2 $ 



E 



{5xj~ 2 



(12) 



where « represents spatial index and n temporal one. 



III. RESULTS 



In order to judge whether the global string network relaxes into the scaling regime, we 
give time development of £, which is defined as 

P = ^/t 2 , (13) 
where \x = 2irrj 2 \n(t/(5£ 1 / 2 )) is the average energy per unit length of global strings. 



A. String identification 



Before obtaining £ and the loop distribution function, we must identify the string seg- 
ment. Since spacetime is discretized in our simulations, a point with $ = corresponding 
to a string core is not necessarily situated at a lattice point. In the worst case, a point with 
$ = lies at the center of a plaquette. Therefore, we use a static cylindrically-symmetric 
solution, which is obtained by solving the equation 

^ + ;f-£-™ = o, (i4) 

with <fr(r,0) = f(r)e ld and the winding number n — 1. The boundary conditions are given 
by 
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/(r) -> |$| min , (r -»• oo) , (15) 
/(0) = 0. (16) 

We require that a lattice is identified with a part of a string core if the potential energy 
density there is larger than that corresponding to the field value of a static cylindrically- 
symmetric solution at r = 5x p h ys /v^2- Then only one lattice within a section of a straight 
string core is identified with a string segment except the case where a point with $ = 
lies at the center of a plaquette. Of course, the real string is not exactly straight but bent 
and more complex. But, as seen in Fig. 0, our identification is good and only one lattice 
within a section of a string core is identified with a part of a string core. In order to solve 
the above equation of motion, we have to use the standard shooting technique, which is the 
repeated process so that it is bothersome to follow the above procedure for each time step. 
Instead, we obtain the solutions every 500 8t and make the fitting formula. For intermediate 
time steps, we judge whether a lattice point belongs to the string segment by comparing the 
potential energy obtained from the simulation and that from the fitting formula. Thus, by 
counting the number of the lattices identified with a part of a string core, we can evaluate 
the total length of strings within the horizon volume, (2t) 3 , from which the energy density 
can be obtained. 



B. Scaling property 

Time development of £ in the cases from (1) to (4) is described in Fig. |3|. We find that 
after some relaxation period, £ becomes a constant irrespective of time with (1)0.99 ± 0.09, 
(2)0.97±0.04, (3)0.90±0.06, and (4)0.90±0.05. They all are consistent within the standard 
deviation. Thus we can conclude that a global string network relaxes into scaling regime in 
the radiation domination. We also show time development of £ in the case (7) in Fig. |4j 
Then £ asymptotically becomes a constant, 0.88 ± 0.07, which is also consistent with the 
above all cases with ( = 10 within the standard deviation. Hence we can also conclude that 
( does not change the essential result. Note that the standard deviation in the case (2) is 
much smaller than the other cases because the box in the case (2) includes more horizon 
volumes at each time. Also, £ seems to oscillate at the early epoch because the homogeneous 
mode of the field oscillates in the radial direction of the potential, which rapidly decays. In 
fact, the period of the oscillation coincides with 2ir times the inverse mass at the potential 
minimum. 

Fig. [| represents the results under the reflective boundary condition, where £ becomes a 
constant irrespective of time with (1)1.77±0.03, (2)1.57±0.04, (3)2.00±0.05, (4)1.30±0.02, 
(5)1.18 ± 0.03, and (6)1.03 ± 0.03. Though all results are consistent with that under the 
periodic boundary condition within the factor two, the former tends to become larger than 
the latter. This is because under the reflective boundary condition, strings are repulsed by 
the boundary and the string near the boundary intercommutes less often than that near 
the center of the simulation box because the partner to intercommute only lies in the inner 
direction of the boundary. The results in large box simulations, (4)- (6), are converging so 
that it is safe to say that if we take the box size larger than 2 3 times the horizon volume, 
the reflective boundary condition has no effect on the results. 
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FIG. 2. Snapshots of a realization at t = 35, 45, 55, 65, 75 in the case (1) under the periodic 
boundary condition. Lattices identified with a part of a string core are shown. 
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FIG. 3. Time development of £ in the cases from (1) to (4) under the periodic boundary 
condition. Symbols(n) represent time development of £. The vertical lines denote a standard 
deviation over different initial conditions. 
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FIG. 4. Time development of £ in the case (7). Symbols(n) represent time development of £. 
The vertical lines denote a standard deviation. 

C. Loop distribution 



We also investigate the loop distribution. Since in our simulations it is judged by the 
potential energy whether a lattice point is a part of the string, the string is not necessarily 
continuous. Therefore, we identify a closed loop as follows; First we select a lattice point 
belonging to the string segment. Then we connect it with the nearest neighbor among lattice 
points belonging to the string segment. We proceed this process one after another until the 
connection returns to the starting lattice point. 

Kibble's one scale model predicts the loop distribution as 



ni(t) 



v 



t*(i + Ktyi 



(17) 



where v is a constant, I is the length of a loop, and the log- dependence of /i is neglected. 
Different from a local string, the dominant energy loss mechanism for global strings is 
radiation of the associated Nambu-Goldstone field ||24j| . We define the radiation power, P, 
as P = Kfi where k is a constant. An example of decay of a closed loop is depicted in Fig. 

s 

We determine whether the loop distributions in the simulation coincides with the above 
predicted loop distribution function. The loop distributions in the case (1) at different times 
(t = 45, 55, 65, 75ij) are described in Fig. [7]. Since long strings are rare, we cut the length 
of loops into bins with the width 5x5x. Also, we divide 300 realizations into 6 groups 
comprised of 50 realizations and we summed the number of loops over 50 realizations for 
each groups. The dot represents the number of loops averaged over 6 groups and the dash 
line represents the standard deviation. They can be simultaneously fitted with the above 



10 




20 40 60 8D 20 40 60 80 

time[t/t.] 



FIG. 5. Time development of £ in the cases under the reflective boundary condition. 
Symbols(n) represent time development of £. The vertical lines denote a standard deviation. 
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formula if one takes k ~ 0.535 and v ~ 0.0865. Fittings for k and v are also given in Fig. 
H Thus, the loop production function as well as the large scale behavior of the string scales 
together for the global string network. 

As an another evidence for the scaling of the loop production, we consider the NG boson 
spectrum radiated from strings. If loops are formed at scales much smaller than the horizon 
distance, there should be significant power of radiated NG bosons for modes corresponding 
to the scale at which loops are formed. For the purpose, we first represent the complex field 
$(£, x) in terms of the radial mode x(t, x ) and the NG boson field a(t, x) as 



V2 



exp 



ia(t, x) 



(18) 



Then the kinetic energy density of NG bosons is given by 

-a(t,xY = - — - 1 — — - 
2 1 ; mt,x)\* 



X 



-Im$(t, #)Re$(t, x) + Re$(t, #)Im$(t, x) 



(19) 



One may wonder if the power spectrum obtained by the Fourier transformation of the 
above kinetic energy density is what we want. But it includes the contribution from NG 
bosons formed at the symmetry breaking as well as that radiated from strings. Also, the 
decomposition of the field into the radial and phase modes is inadequate for the lattice point 
near the string segment. 

Then, we evaluate the average energy density of NG bosons radiated in the period be- 
tween t\ and t 2 , p[ti,t 2 ). For the purpose, we subtract the redshifted kinetic energy at ti 
from the kinetic energy at t 2 since emitted NG bosons damp like radiation. Thus, p[ti,t2] 
is given by 



1 r 

p[ti,h] = y J d 3 xp[t 1 ,t 2 ;x] 



v ■ dx 

1 r d 3 k 



-a(t 2 ,x) - -a(h,x, 



VJ (2vr) 3 



- \a k {t 2 )\ 2 - - 1^^; 



■f \ 2 

|2 / J l 



d 3 k _ . 



(2tt 



dk 
2^2 



where V is the simulation volume and a(t, x) = J ■^^<^h(t) exp(ifc • a;). 



2J , 



(20) 



Furthermore, to avoid contamination of string cores to the spectrum of emitted axions, 
we divide the simulation box into 8 cells and stock the field data of a cell if there are no string 
cores in that cell between t\ and t 2 - Over all such cells, we average power spectra of kinetic 
energy of axions obtained through Fourier transformation. We follow the above procedure 
between t\ = 65ti and t 2 = 75tj for the case (3) under the periodic boundary condition, which 
is the highest resolution simulation, but with the zero temperature potential V e s[^,T = 0] 
after t = 20tj because the decomposition of the field becomes well-defined. One may suspect 
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FIG. 7. Loop distributions at t = 45, 55, 65, 75 are depicted. The number is summed over the 
box size (128(5x) 3 ) and 50 realizations for each groups. Bins are cut every 5x5x. 
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FIG. 8. Fittings for n and v. The cross point represents the best fit values for k and v. The 
solid circle denotes 68% C.L. 



that this spectrum is dominated by the kinetic energy of NG bosons associated with global 
strings (string NG bosons) rather than that of free NG bosons radiated from global strings. 
But this is incorrect for the following reason; The total energy of string NG bosons is almost 
as much as that of free NG bosons. However, the contribution to the energy of string NG 
bosons is dominated by the gradient energy because the kinetic energy of string NG bosons 
is much smaller that their gradient energy by the factor v 2 (v is the velocity of the string 
core and i;<l except just before the disappearance of the loop). On the other hand, the 
gradient and the kinetic energy of free NG bosons is the same. Therefore the kinetic energy 
of free NG bosons is much larger than that of string NG bosons, that is, our spectrum 
is dominated by the former. Also, even if it contributed, the kinetic energy of string NG 
bosons would decay in proportion to t~ 2 , so that its contribution has been removed from 
the spectrum in our method as done in Eq. 



The result is depicted in Fig. ^j, which has already been given in |23| in the different 
context. The spectrum of emitted NG bosons is highly peaked at the horizon scale. Thus, 
loops are formed not at much smaller scale than the horizon distance but around the horizon 
scale. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this paper we gave a comprehensive investigation on the evolution of the global string 
network in the radiation dominated universe by use of the numerical simulation based on 
the complex scalar field model which spontaneously breaks the U(l) symmetry. 

In order to decide whether the large scale behavior of the global string network goes 
into the scaling regime, we followed time development of the scaling parameter £ which 
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Spectrum of NG bosons emitted between t = 65t. and t = 75L 
,0.01 1 1 1 1 1 1 — ^ 1 



LO 
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wave number L(k/5k)~ LkJ 

2tt 



FIG. 9. Filled dots represent the power spectrum, p^, which is already averaged over the 
direction of k and multiplied by the phase-space factor as defined in Eq. (]20j) . Bins are cut every 
55k. k = 645k corresponds to string cores. 



characterizes the average energy density of global strings. We found that £ is almost constant 
irrespective of cosmic time under both the periodic and the reflective boundary condition. 
All the results are consistent within the factor two. But £ under the reflective boundary 
condition tend to become larger than that under the periodic boundary condition. This can 
be understood as follows; Under the periodic boundary condition, there are no infinite strings 
and strings with two boundary points on the opposite planes always intercommute with the 
partner so that £ tends to become small. On the other hand, under the reflective boundary 
condition, strings are repulsed by the boundary. Furthermore the string near the boundary 
intercommutes less often than that near the center of the simulation box because the partner 
to intercommute only lies in the inner direction of the boundary. Thus, £ tends to become 
large. Therefore, £ in the real world should lie between that under the periodic boundary 
condition and that under the reflective boundary condition. Considering all the cases under 
the periodic boundary condition and the cases (4)-(6) of the large box simulations under the 
reflective boundary condition, it is safe to say that £ ~ (0.9 — 1.3). 

You should note that £ is much smaller than that of a local string, which is of the order 
of ten. This is mainly because global strings can intercommute more often than local strings 
since an attractive force proportional to the inverse separation works between a global string 
and a global anti-string. 

We have also investigated the loop distribution. It can be well fitted to that predicted by 
the one scale model if we take v ~ 0.0865 and k ~ 0.535. Thus, the loop production grows 
with the horizon distance and we did not observe the small scale structure. This is because 
the Nambu-Goldstone(NG) boson radiation from strings is so efficient that the small scale 
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structures on strings are damped out. The damping scale is typically Kt, which is near the 
horizon distance in contrast with Gfit <C t (G : the gravitational constant) if at all for the 
local string network. Bennett |J showed that unless produced loops with the length of the 
horizon distance often self-intersect and fragment into smaller loops with the typical length 
smaller than the horizon distance, the reconnection rate is large enough to prevent scaling. 
In the global string network, parent loops do not fragment into smaller loops with the typical 
length smaller than the horizon distance, instead they rapidly shrink through radiating NG 
bosons so that the reconnection rate of parent large loops becomes small and scaling should 
be kept. 
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